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Abstract. Recently several quantum search algorithms based on quantum walks 
were proposed. Those algorithms differ from Grover's algorithm in many as- 
pects. The goal is to find a marked vertex in a graph faster than classical algo- 
rithms. Since the implementation of those new algorithms in quantum computers 
or in other quantum devices is error-prone, it is important to analyze their ro- 
bustness under decoherence. In this work we analyze the impact of decoherence 
on quantum search algorithms implemented on two-dimensional grids and on 
hypercubes. 

Resumo. Recentemente, vdrios algoritmos qudnticos de busca baseados em 
passeios aleatdrios qudnticos foram apresentados. Estes algoritmos sao difer- 
entes do algoritmo de Grover em vdrios aspectos. O objetivo e encontrar um 
vertice marcado em um grafo mais rdpido do que algoritmos cldssicos. Uma vez 
que a implementagao destes novos algoritmos em computadores qudnticos ou 
qualquer outro dispositivo qudntico estd sujeito a erros, e importante analisar a 
robustez em relagao a descoerencia. Neste trabalho, analisamos a descoerencia 
dos algoritmos qudnticos de busca em malhas bi-dimensionais e hipercubos. 

1. Introduction 



After Grover's seminal paper JGrover (1996)H , it is known that a quantum computer can 
search an element in a non- structured database quadratically faster compared to a clas- 
sical computer. Recently many other search algorithms were developed based on the 
discrete-time quantum walk model [ Kempe (2003)[ . They differ from Grover's algorithm 
in many aspects and may be better suited for practical implementation ]SKW (2003 )[ 



AKR (2005)|. 



Quantum walks can be seen as a quantum counterpart of classical random walks 
[ |Kempe (2003)[ |. They have been successfully applied in several quantum algorithms 
QAmbainis (2004)j |Szegedy (2004)1 |Magniez, Santha and Szegedy (2007)). It is possi- 



ble to describe an abstract search algorithm JAKR (2005)[ in a generic regular graph 
by using a discrete-time quantum walk with a modified coin operator to search for a 



marked vertex. Grover's algorithm can be seen as the simplest example of this proce- 
dure. The Shenvi-Kempe-Whaley (SKW) algorithm can be seen as an implementation 



of the abstract search algorithm to an n— dimensional hypercube QSKW (2003)| and the 
Ambainis-Kempe-Rivosh (AKR) algorithm is an application to two-dimensional grids 
QAKR (2005)1 . 



Decoherence and gate imperfections are unavoidable side-effects in any imple- 
mentation of a quantum computer. One approach to deal with this problem, with consid- 
erable overhead in quantum resources, consists in using redundant encoding and several 
layers of error-correction codes. Another approach consists in designing algorithms that 
are robust against certain types of errors which may be dominant in a given implementa- 
tion. This requires a detailed knowledge of the effect that different kinds of noise have on 
the performance of the algorithm. It seems likely that a real quantum computer will take 
advantage of both approaches. Decoherence in quantum walks has been considered in 
previous works, see for instance [ |Kendon and Tregenna (2003)[|Romanelli, et. al. (2004)[ 



Alagic and Russell (2005) |. 



In this work, we study the effect of three different noise models on the AKR and 
SKW algorithms without quantum error correction. In two of them we assume phase 
errors affecting the modified coin operator. In the other case, we assume broken-link 
imperfections affecting the shift operator. The influence of errors in the AKR algorithm 
has not been analyzed before. This work presents the first results in this direction. Li and 



collaborators [ |Li, Ma, and Zhou (2006)] have analyzed the effect on the SKW algorithm 
of gate errors in the coin operator assuming a perfect shift operator. In this work we 
complement those results by analyzing errors in the shift operator. 

The structure of the paper is as follows. In Section 2, we review Grover's algo- 
rithm and introduce the notation used in this work. In Section 3 we review some basic 
concepts of search algorithms based on discrete-time quantum walks. In Section 4, we 
describe the noise models studied in this work and review the broken-link type of noise 



[Marquezino, Portugal, Abal and Donangelo (2008) |. In Section 5, we review the SKW 



Algorithm. In Section 6 we review the AKR Algorithm. In Section 7, we describe the 
methodology used to study each of the noise models and give the results of the numerical 
simulations we performed for the SKW algorithm. In Section 8, we present the results of 
the numerical simulations for the AKR algorithm. In Section 9, we present our conclu- 
sions. 

2. Grover's Algorithm 



In this section we briefly review Grover's algorithm JGrover (1996)] and the notation used 



in this paper. For further details see [Nielsen and Chuang (2000) |. Consider a quantum 



computer of n quantum bits (qubits). Quantum Mechanics tells us that its state is de- 
scribed by a unit vector in a vector space of dimension N = 2 n . The simplest orfhonormal 
basis for this vector space is {|0), |1), • • • , \N — 1)}, where \i) is a vector of N entries all 
of them zero except the entry i + 1 which is 1. This basis is called computational basis. 
Suppose that the state of the quantum computer in a given instant of time is the vector 
l^), then 

N-l 
i=0 



where the coefficients ai{ are complex numbers that must obey the constraint 



E h 2 = L 



(2) 



i=0 



It is possible to prepare the quantum computer at the beginning of an algorithm in any 
state The algorithm must be a sequence of applications (multiplications) of N x N 
unitary operators (matrices) U\, ■ ■ • ,U t . So, at time t, the state of the quantum computer 
is 



A matrix U is unitary if UW = I, where W is the transpose conjugate of U and I is the 
identity matrix. 

Since the quantum computer is a physical system, one can perform measurements 
to determine the state Quantum Mechanics tells us that the result of a complete 
measurement in the computational basis will be i with probability | otj | 2 . The measurement 
does not allow one to find the coefficients a*, that are needed to describe the vector \ip). 
Instead, one gets a random number in the set {0, • • • ,N — 1} with the corresponding 
probability distribution {|« | 2 , • • • , |aAr_i| 2 }. A complete measurement means that all 
qubits are measured yielding either or 1 each. One gets a binary number that is converted 
to the decimal notation. A partial measurement consists in measuring a fraction of the 
qubits. 

A search algorithm, such as Grover's algorithm, aims to determine whether an 
element i belongs to non-structured database or to determine the position of an element 
in a non-sorted database. The easiest way to pose this problem is in the following form. 
Suppose that the domain of function f is {0, ■ ■ ■ , N — 1} and the image is 



Suppose that we ask a friend to implement the function / in a classical computer and 
there is only one value of x such that f(x) = 1. We can obtain the image of any input 
value by employing /. What is the complexity of the best algorithm that finds the value 
of x such that f(x) = 1? The complexity in this case is measured by the number of times 
we employ the function /. If we do not know any equation for / nor any details of the 
implementation of /, the only way to find out that f(ia) = 1 is through an exhaustive 
search. The complexity of the best classical algorithm is O(N). 

In a quantum computer, the implementation of function / must be performed 
through a unitary matrix, which we call Uf. The definition of Uf is 



\i;) = U t ---U 1 \i). 



(3) 




(4) 




(5) 



We suppose that the quantum computer has n + 1 qubits. So the output of the function / is 
added to a second register of 1 qubit. We again ask our friend to implement U f, this time 
in a quantum computer. What is the complexity of the best quantum algorithm that finds 



the value of x = io? Grover's algorithm requires O(yN) applications of Uf to determine 
the value io with a very small margin of error. 

The idea behind the algorithm is to start the quantum computer in a known state, 
then apply a sequence of unitary matrices that results in a state that has a large overlap 
with state |z'o). If the state of the quantum computer has a large overlap with \i ), that 
is (ip\io) ~ 1, the result of the measurement will be io with high probability. The notation 
(ip\ means the transpose conjugate of \ip). 

Grover's algorithm works as follows. Prepare a (n + l)-qubit quantum computer 
in state \ip ) = |0)|1). Apply #®( n+1 ) to |^ ). H is the Hadamard matrix and <g> is 



the tensor product, see [Nielsen and Chuang (2000) |. The result is \ip\) = \s)\—), where 



-) = (|0)-|l))/ v / 2and 

N-l 







Now apply t/TIv^l t0 where 

U=((2\s)(s\-I N )®I 2 )U f . 
In the last step, measure the first register to obtain i with probability 1 — 0(1 /N). 

3. Quantum Walk based Search Algorithms 

Quantum walks generalize the concept of classical random walk. A classical walk is a 
prescription of how to move, conditioned to the value of a random variable. If the walker 
lives in a regular graph of degree d, the random variable must have d values, usually with 
the same probability l/d(a balanced walk). The edges of the graph incident to a vertex v 
must have labels from to d — 1. If the walker is in vertex v and the result of the random 
variable is j, than the walker moves to the vertex v' that is connected to v by an edge of 
label j. This procedure is repeated again and again. The result is a random walk on the 
graph. In a one-dimensional lattice, one can toss a coin and move to the right if the result 
is heads, or to the left if the result is tails. 

In a quantum setting, both the toss of a coin and the shift of the walker must 
be performed by unitary operators. In a regular graph of degree d, the vector space 
where the walk takes place is He <8> Hy, where He is the Hilbert space spanned by 
{|0), • ■ ■ ,\d — 1)} representing the coin space and Hy is the Hilbert space spanned by 
{|0), • ■ ■ , | V — 1)} representing the vertex space, where V is the number of vertices. The 
usual form of the evolution operator is 

U = S(C® I). (7) 

Here, C is a d x d matrix that acts only on the coin space, I is the identity in the vertex 
space and S is the shift operator given by 

s\m = \j)w, (8) 



where v' is the vertex that is connected to v by edge j. Note that the coin operator C is 
the same for all vertices. The walker starts at some initial configuration called \ipo) and at 
time t its state is U*]^}. 



Quantum walk search algorithms are based on a modification of some standard 
quantum walk given by U. The analysis of the new walk depends on what happens in 
the original non-modified walk. Suppose that we would like to search for vertex vq. That 
vertex must be marked somehow. One marks this vertex using a modified coin operator. 
The new coin operator must distinguish the marked vertex from the rest The new coin is 
defined by 

C = (-/) ® \v )(v \ + C <g> (/ - \v )(v \). (9) 

This new coin operation C defines a new evolution operator given by U' — S C. 

The most used coin is called Grover's coin and it is the real unitary operator far- 
thest from the identity [ |Moore and Russell (2002)1 . It is defined as C = 2\s)(s\ — I, 



where \s) is the uniform superposition, 

For this coin, all directions have the same weight. Using Grover's coin in Eq. ©, one 
obtains 

U' = U(I-2\s,v )(s,v \). (11) 
It is possible to perform a very general analysis of search algorithms on graphs if one 



demands some properties from U [AKR (2005) | . These properties are: (1) U must be 



a real unitary matrix, (2) U must have only one eigenvector with eigenvalue 1, and (3) 
this eigenvector must be the initial state of the algorithm. This abstract search algorithm 
works as follows. Suppose that U\4>q) = |0 O ) an d let exp(za) be the eigenvalue of U' with 
the smallest angle among all eigenvalues of U' . If U' has eigenvalue 1, the initial condition 
and the evolution of the walk must be in a space orthogonal to the eigenspace associated 
with eigenvalue 1. The algorithm consists in applying U' ^1 to |0o) an d measuring the 
vertex space. 

Grover's algorithm is the simplest example of form (fTTT) . Two new search al- 
gorithms in this framework have been analyzed in detail. The first one is a search on 
hypercubes (SKW), which we describe in Sec. [5J The second one is a search on two 
dimensional grids (AKR), which we describe in Sec. [61 



4. Decoherence models 

In an actual physical implementation, operators are error-prone. It is important to deter- 
mine the robustness of an algorithm to errors in its implementation. In quantum search 
algorithms, there are three key operators: the original coin (C), the coin used in the 
marked vertex (— /) and the shift (S). We analyze the impact that errors on each of these 
operators have on the algorithm's performance. 

Phase errors on the coin operator affecting the marked node can be implemented 
by replacing — J by 

C V0 (9) =e^ +e h (12) 

with 9 E [— 7r, 7r], The perfect coin operator — / on the marked vertex is recovered for 
8 = 0. We say that the error is systematic when the phase error 9 is constant in each step 



(model I), and that the error is random when the phase error 6 in each step is a Gaussian 
random variable with zero mean and standard deviation a (model II). 



Phase errors on the coin operator for the unmarked nodes can be implemented by 
rewriting C as 

C{6) = I-(l-e i{7T+e) )\s)(s\ (13) 
with 9 E [—tv, 7r]. The Grover coin operator is recovered for = 0. 

The effect of phase errors in the original Graver's algorithm was analyzed by 
Long and collaborators and, later on by Shenvi and collaborators [ |Long, et. al. (2000)} 



Shenvi, Brown and Whaley (2003) |. The latter authors investigated the importance of the 
scaling of phase errors with the size of the database, N . In a recent work, Li and col- 
laborators [ |Li, Ma, and Zhou (2006)| studied the effect of an imperfect C on the SKW 



algorithm. The operators —I (acting on the marked vertex) and S were supposed to have 
no errors. 

Errors in the shift operator S can be implemented by randomly opening links 
between connected vertices with probability p per unit time (model III). This broken- 
link noise model has been previously considered for a quantum walk on a line and on 
a plane [ |Romanelli, et. al. (2004)[ |01iveira, Portugal and Donangelo (200~6)fl , and on the 



hypercube [Marquezino, Portugal, Abal and Donangelo (2008) |. To implement this kind 



of error we generalize the shift operator S such that no probability flux is transferred 
across a broken link. This modified shift operator is unitary for any number of broken 
links in the lattice. At each time step, the topology of the graph is defined, opening each 
link with probability p and performing the shift to the neighboring vertex only if the link 
is not broken. The original S operator is recovered for p = 0. 

5. SKW Algorithm 

The quantum search on n-dimensional hypercubes has a Hilbert space He ®Hp, where 
He is a n-dimensional Hilbert space associated with a "quantum coin" and Hp is a 2 n - 
dimensional Hilbert space associated with the vertices of the hypercube. A basis for He 
is {\d)}, for < d < n — 1 and a basis for Hp is {|x n _ix„_ 2 • • • x )}, for binary x d . In 
a hypercube two vertices are connected if, and only if, the corresponding binary strings 
differ by one bit. 

The generic state of the quantum walker is 

n-1 2" 

i*w) = EEfeWi d ' f )' (14) 

d=0 x=0 

where ijjd;x(t) £ C and Yld=o l^tyzWI 2 = 1- The evolution operator for one step of 
the walk is U — S (C (g) I), where S is the shift operator and C the coin operator acting 
in He and I is the identity in Hp. The shift operator can be written as 

n-l 

S = \d, x ® ed){d, x\, (15) 

d=0 x 

where © is the bitwise binary sum between binary vectors, and is a null vector except 
for a single 1 entry in the <i-th component. The coin operator is given by 

C = 2\s c )(s c \ -I, (16) 



where \s°) = Yld=o M) * s the uniform superposition over all directions. 

In the SKW search algorithm, the initial state is the uniform superposition of coin 
and spatial eigenstates, i.e., |^ ) = \ s °) ® l sP )' where \s p ) is the uniform superposition 
over the 2" nodes of the hypercube. Note that a single searched state is contained in \s p ) 
with an amplitude l/\/2™- The initial state |^o) is an eigenstate of U with eigenvalue 1. 



The quantum search (SKW), as proposed by ]SKW (2003)1 , is based on a modified coin 
operator of the form 

C' = C®/-(/ + C)<g)|0)<0|. (17) 

Without loss of generality, we assume that the searched node is labeled |0). The detailed 
analysis of this algorithm shows that after |V2"^ = 0(\/N) iterations, a measurement 
of the position of the walker yields the marked state with success probability \ — 0(l/n) 



[SKW (2003)| 



6. AKR Algorithm 

The quantum search on a \fN x \f~N grid has a four-dimensional Hilbert space He for 
the coin and a ^-dimensional Hilbert space Hp for the N lattice sites. A basis for He 
is {\d,j)}, for < d < 1 and < j < 1. Variable d sets the direction of the walk, 
d = for horizontal shift and d = 1 for vertical shift. Variable j sets whether the 
walker moves forwards (j = 0) or backwards (j = 1). A basis for Hp is {|no, ^i)}, for 
< n 0l ni < y/N. The boundary conditions are periodic. 

The generic state of the quantum walker is 

1 

Yl ^d^nMdJ^nt). (18) 

d,j=0 no,ni=0 

The action of the shift operator on the computational basis is 

S\d,j)\n , ni ) = \dJ®l)\no + (-l) j 5 dfl , ni + (-1)^}, (19) 

where © is the binary sum. Notice that there is an inversion from backwards to forwards 
and vice-versa after the action of the shift operator. This is a modification of the standard 
shift operator and it is absolutely necessary to have a quadratic speedup over the optimal 
classical algorithm. 

The initial state is the uniform superposition of coin and spatial eigenstates, i.e., 
\ijj ) = \s c ) ® |s p ), where \s p ) is the uniform superposition over the N sites of the grid 
and 

d,j=0 

The coin operator on unmarked vertices is given by Eq. (fT6l) and the modified coin oper- 
ator is 

C' = C®I- (I + C) ® |n ,ni)(n ,ni|, (21) 

where (n , ni) is the marked vertex. The evolution operator is U' = S C and must be 
applied 0(\JN\og N) times. The overlap between the final state and the marked ver- 
tex is 0(l/v / loglV). In order to improve the probability of finding the marked vertex, 



it is necessary 0(y/\ogN) rounds of the algorithm yielding an overall complexity of 
0(viVlog N). A very recent paper [ |Tulsi (2008)H improved the overall complexity to 

0(^N log N). 



7. Results for SKW Algorithm 

Fig. Q] shows the probability of finding the walker at the marked node as function of the 
number of steps, for each noise model. In the left panel we compare the results for the 
ideal case, without noise, with those for noise models I and II. In the right panel we 
compare the ideal case with the evolution under noise model III. All plots correspond 
to a hypercube of dimension n = 8. For model I, we took a phase error 9 = 0.3, the 
standard deviation in model II was a = 0.3 and the probability of broken links per unit 
time (model III) was p = 0.02. Note that the peak of probability for the systematic 
error (model I) occurs earlier than in the ideal case with zero noise. The behavior of 
the algorithm under noise from model I affecting the coin operator at the marked node 
(—1) is similar to the one observed in [ |Li, Ma, and Zhou (2006)] , in which the operator 
C was affected. Although random models II and III correspond to different kinds of noise, 
they result in similar patterns. In both cases the first maximum in the probability occurs 
approximately at the same number of steps as in the case with no noise, and it reaches a 
lower value. Subsequent peaks undergo a gradual attenuation with the number of steps s. 
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Figure 1. Left panel: probability at the marked vertex as a function of the number 
of steps s comparing the ideal case with both systematic [6 = 0.3) and random 
(<t = 0.3) errors. Right panel: the same for broken-link errors with p = 0.02. 



The stopping time in the case without noise corresponds to the first maximum and 
this point is ^V2 n ~ 1 ~ 18, for n = 8. In the presence of noise it is better to stop the 
algorithm before this point and to rerun to find the correct result. If the probability to 
obtain the correct result in one run is p, then the expected number of times of trials is 
1/p. If the computational complexity of one run is O(VN) then the overall complexity is 
0(\fN /p). If p does no depend on N, it is not going the change the complexity. Let us 
define the algorithmic cost c(s) as the overall cumulative number of steps needed to find 
the desired state, 

c(s) = -, (22) 

Ps 

where s the number of steps before the final measurement is taken in one run of the 
algorithm. In Fig.[2]we show the cost function, c(s), for the different noise models. In the 
case of a systematic phase error, the cost function has a well-defined minimum at s fa 10. 
It is clearly convenient to stop the algorithm before the peak probability is reached, in 



either the ideal case or in the case with noise. For the other noise models, and also in the 
case without noise, the cost function has a very shallow minimum, and it increases very 
slowly with step number after its minimum. However, even in these cases, these results 
suggest that it is advantageous to stop the algorithm before the noiseless probability peak 
is reached, and to repeat it more times, as needed. 
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Figure 2. Cost c(s) from Eq. (|22), vs. number of steps, for the noiseless search 
algorithm and for the algorithm with the three different noise models described 
in the text. The hypercube considered has dimension n = 8. 

In Fig. [3] we observe the probability of reaching the marked vertex as a function 
of the noise strength, as represented by the three upper curves. In the left (right) panel 
we have the results for model I (model II). Note that the optimal phase is 9 = or 
(7 — 0, i.e., when —I is used in the marked node as in the standard algorithm. The 
plot also shows that the algorithm is very sensitive to noise from operational errors, if 
no error correction code is used. The three lower curves represent the highest probability 
among the unmarked vertices. We observe that, as the phase error increases, the difference 
between the maximum probabilities at marked and unmarked nodes becomes smaller. In 
this case, we cannot distinguish the right solution and the algorithm is no longer useful. 
The noise generated by systematic errors (left panel) seems to play a more significant role 
in the algorithm than the noise generated by random errors (right panel). 
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Figure 3. Left panel: Results for model I. Right panel: Results for model II. Up- 
per three curves: Highest probability at marked vertex as a function of the noise 
strength parameter for three values of the dimension of the hypercube n. Lower 
three curves: Highest probability at the unmarked vertices, using the same con- 
vention for the dependence of the dimension of the hypercube. 

In Fig. H] we have the results for model III. In the left panel, we observe the highest 
probability at marked vertex as a function of the broken-link rate p, as represented by the 
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three upper curves. The three lower curves represent the highest probability among the 
unmarked vertices. In the right panel, we have the highest probability at the marked vertex 
as a function of the dimension n of the hypercube. In this case, the probabilities decay 
as the dimension of the hypercube is increased, which is similar to the result obtained in 
[|Li, Ma, and Zhou (2006)] for noise affecting the coin operator at unmarked sites. 
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Figure 4. Results for model III. Left panel: similar to Fig. [3] as a function of the 
broken-link rate. Right panel: highest probability at marked vertex as a function 
of the dimension n of the hypercube. 



In order to estimate how errors of model II change the complexity of the algorithm 
we use a formula that scales with N in the form 9 = 1/N 5 in Eq. (PT21) . In Fig. \5\ we 
plot the scaled cost, which is the logarithm to base N of the algorithmic cost given by 
Eq. (1221) . against error parameter 5 for several values of n. Recall that the complexity 
of SKW algorithm is 0(N 0,5 ) and its success probability is 1/2 — 0(l/n). Hence, for 
large values of 5 and N, we should obtain a scaled cost close to 0.5, corresponding to 
the complexity of the noiseless SKW algorithm. Our plot shows a scaled cost close to 
0.6, which is consistent with the values of iV considered. This means that for 6 > 1, the 
SKW algorithm with error has the same complexity of the noiseless SKW algorithm. For 
5 < 1, the noise rate increases and the algorithm gradually looses efficiency in relation 
to the noiseless search. For 5 ~ —0.1, the scaled cost is close to 1, which means that the 
quantum algorithm has the same complexity of the classical search, O(N). For S < —0.1, 
the scaled cost is higher than 1, which means that the performance of the quantum search 
is worse than the classical search. 
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Figure 5. Logarithm (base N) of the algorithmic cost as a function of scaling 
parameter S for model II comparing different dimensions. 



8. Results for AKR Algorithm 



The behaviour of the maximum probability at the marked node in the AKR algorithm 
follows a similar pattern as in the SKW algorithm. The main difference being that for 
AKR the maximum probability decreases as N increases, while for SKW the maximum 
probability remains close to 1/2. The numerical results for the cost in the AKR algorithm 
also show that in the presence of imperfections it is better to stop the algorithm before the 
theoretical stopping time. 

Fig. [6] shows the maximum probability at the marked node as function of the noise 
strength. In the left (right) panel we have the results for model I (model II). This figure 
should be compared to Fig. |3] Note that the number of nodes of the grids corresponds 
to the number of nodes of the hypercubes. The curves for the AKR algorithm are very 
similar to those for the SKW algorithm and we draw similar conclusions for both cases. 
The main difference is the distance between the curves, a consequence of the fact that in 
the AKR algorithm the maximum probability at the marked node drops when we increase 
N. 




Figure 6. Left panel: Results for model I. Right panel: Results for model II. High- 
est probability at marked vertex as a function of the noise strength parameter for 
three values of the dimension of the grid. 



In Fig. U\ we show the results for model III. In the left panel, we observe the 
highest probability at the marked vertex as a function of the broken-link rate p. In the 
right panel, the horizontal axis is in log scale. These results should be compared with 
Fig. IU The probability drops faster in the AKR algorithm than in the SKR. This was 
predicted in JAKR (2005)[ , where it is shown that the probability at the marked node 
scales as 0(1/Vlog N). 

In Fig. [8] we plot the scaled cost log N c(s) against 5, which is the analogue of 
Fig.[5]for the AKR algorithm. Recall that the cost in the AKR algorithm is 0( A 5 log N) . 
Hence, for large values of 5 and N, we should obtain a scaled cost a little bit above 0.5, 
corresponding to the complexity of the noiseless AKR algorithm. The scaled cost is not 
exactly the power of N because the cost has the term log N. Our plot shows a scaled cost 
close to 0.8, which is consistent with the values of N considered. From the figure we see 
that for 5 > 1/2, the AKR algorithm with error has the same complexity as the noiseless 
AKR algorithm. For 5 < 1/2, the noise rate increases fast enough such that the algorithm 
looses efficiency in relation to the noiseless search. When we decrease 5, the scaled cost 
approaches 1 , which means that the quantum algorithm has the same complexity of the 
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Figure 7. Results for model III. Left panel: highest probability as a function of the 
broken-link rate. Right panel: highest probability at marked vertex as a function 
of the dimension log TV of the \fN x Vn grid. 



classical search, O(N). For 5 < 0, the scaled cost is higher than 1, which means that the 
performance of the quantum search is worse than the classical search. Note that 5 = 1/2 
is the transition point in the AKR algorithm, while 5 = 1 is the transition point in the 



SKW algorithm. For comparison, note that |Shenvi, Brown and Whaley (2003) | obtained 



5 = 1/4 as the transition point in the original Graver's algorithm. 
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Figure 8. Logarithm (base N) of the algorithmic cost as a function of scaling 
parameter S for model II comparing different dimensions in AKR algorithm. 



9. Conclusions 

In this work, the effect of imperfect quantum operators on search algorithms based on 
quantum walks has been investigated. We have considered both systematic and random 
phase errors in the coin operator. The effect of randomly broken links affecting the shift 
operator has also been considered. This kind of error directly affects the spatial propa- 
gation of the walker. We have considered the search of a marked vertex on hypercubes 
(SKW algorithm) and two-dimensional grids (AKR algorithm). 

For the SKW algorithm, we found that the overall effect of noise in the coin op- 
erator for the marked node is similar to that for the imperfect coin operator acting on 
unmarked nodes considered in QLi, Ma, and Zhou (2006)1 - There are also many similari- 



ties with the AKR algorithm. The overall qualitative effect of noise seems to be similar in 
all search algorithms considered. On the other hand, we have obtained quantitative results 
for the tolerance of the algorithms to errors. 



In the context of Grover's search algorithm, it was shown analytically that phase 
errors 9 scaling as 1 /N 5 , for 5 < 1/4, modify the complexity of the algorithm to O^N 1 ^ 25 ) 



[Shenvi, Brown and Whaley (2003) |. For 5 > 1/4, the complexity of Grover's algorithm 
with errors is equal to the complexity of the noiseless case, 0(\/N). If 5 < 1/4, the 
advantage over a classical search, O(N), is progressively reduced. If 5 = 0, i.e., constant 
error, the complexity of Grover's algorithm is equal to the complexity of the classical 
search. We have also considered this setup in the context of the SKW and AKR search 
algorithms. Our numerical simulations show that the transition point for AKR is around 
5=1/2 and for SKW it is around 5 = 1. For 5 below these threshold values, the algo- 
rithm gradually looses efficiency until becoming worse than the classical case at around 
5 = 0. 

Our numerical results show that it is possible to improve the efficiency in all cases 
(with or without noise) if we stop the search algorithm before the number of steps pre- 
dicted theoretically. In this case, more than one round of the algorithm is needed, keeping 
the overall cost smaller than when using the theoretical stopping point. 
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